High-bias stability of monatomic chains 
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For the metals Au, Ft and Ir it is possible to form freely suspended monatomic chains between 
bulk electrodes. The atomic chains sustain very large current densities, but finally fail at high bias. 
We investigate the breaking mechanism, that involves current-induced heating of the atomic wires 
and electromigration forces. We find good agreement of the observations for Au based on models 
due to Todorov and coworkers. The high-bias breaking of atomic chains for Ft can also be described 
by the models, although here the parameters have not been obtained independently. In the limit of 
long chains the breaking voltage decreases inversely proportional to the length. 
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I. INTRODUCTION 

As the down-scaling of electronic components is ad- 
vancing at an impressive rate, the general interest in dis- 
sipation and thermal transport at comparable scales is 
increasing. Technical applications using these small scale 
devices will have higher and higher densities of compo- 
nents making heat generation one of the biggest problems 
to overcome. We, therefore, need to develop a better 
fundamental understanding of the properties of matter 
at small size scales. A recent review^ on thermal trans- 
port at the nanoscale shows how fast this field has been 
growing during the last few years. It also shows that at 
scales of the order of several nanometers the quantum 
behavior of matter starts to have a major impact on the 
thermal behavior of the sample. As an example, Yao et 
al^ have shown that the current through a single-wall 
carbon nanotube is saturating as the voltage reaches val- 
ues of several volts, a behavior which is not expected 
from a simple down-scaling of classical effects. 

Here we study the thermal behavior of metals at the 
smallest length-scales, below 1 nm. A review of results in 
this field has recently appeared^. In recent experiments it 
was shown that some metals have the property of sponta- 
neously forming monatomic chains just before rupture^. 
Although atomic chains are intrinsically metastable, ex- 
periments at high bias voltages showed they have a sur- 
prising current-carrying capacity^, and current densities 
up to 10^ A/mm^ could be reached. Currently the 
lengths of these chains is limited to about seven atoms. A 
better understanding of dissipation and relaxation phe- 
nomena may lead the way to fabricate longer chains. 
With these monatomic chains being nearly ideal one- 
dimensional conductors this would provide means to fur- 
ther study the fundamental properties of conductors in 
one dimension. 

In the experiments described below we probe the sta- 
bility of many individual monatomic chains by looking 
at their rate of survival at high bias voltages. Although 
these nanowires due to their small size are ballistic, high 



bias results in a small probability for inelastic scatter- 
ing of the electrons to phonons. Energy dissipation has 
been demonstrated to come into play when eV for the 
bias voltages applied is raised above the characteristic 
zone-boundary phonon energy^. Energy dissipation will 
take place on the scale of the electron-phonon scattering 
length (typically ^lOOnm), but the number of scatter- 
ing events per volume will be largest inside the atomic 
wire itself due to the large current density. As a conse- 
quence the average lattice vibration energy at the narrow- 
est cross-section will be larger than anywhere else in the 
contact, which can be viewed as an increased local effec- 
tive temperature. Toff. This elevated local temperature 
can lead to the rupture of the contact and is therefore an 
important parameter. 

There have been several experimental attempts in mea- 
suring this local effective temperature inside metallic 
point contacts, but it has been difficult to obtain quan- 
titative results. Van den Brom et alJ have exploited 
the switching rate between two different configurations 
in atomic sized contacts. When the energy barrier be- 
tween the two configurations is low enough (which is of- 
ten found on the verge of an atomic reconstruction) the 
local temperature can activate the jumps between the 
two levels, and the switching rate is determined by the 
ratio of the effective barrier and the thermal energy. To 
obtain the effective temperature from this switching rate 
an estimate of this barrier is required, which is difficult in 
view of the unknown atomic configurations. In order to 
circumvent this problem the switching rate between the 
two states was followed while changing the bias voltage 
and the bath temperature independently*. A comparison 
between the changes in the switching rates in these two 
cases would provide a calibration of the relation between 
the bias voltage and the effective temperature. How- 
ever, the thermal expansion in the electrodes due to the 
changing temperature changes the atomic configuration 
considerably, making it difficult to draw definite conclu- 
sions. 

Here we discuss measurements of the fracture rate of 
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atomic chains as a function of the voltage bias. The ad- 
vantage over previous work on T^g in atomic-sized con- 
ductors is that they are free of lattice imperfections and 
impurities and the structure is known, except for the con- 
figuration at the end points. The energy barrier that 
has to be considered in the process is therefore less com- 
plex. Moreover, it is possible to consider chains of vari- 
ous lengths, which provides us with a control parameter 
of the point contact geometry. 

The paper is organized as follows. After some details 
on the experimental techniques we show in section IIIII 
that the breaking of an ensemble of chains cannot be 
described by a simple activated process. We argue (sec- 
tion ^3 that there are a number of parameters involved 
that lead to a distribution of Boltzmann factors and we 
show that the observed distribution in breaking voltages 
can be fairly well described without adjustable param- 
eters. One of the assumptions, here, is that there is a 
variation in effective barriers to breaking as a result of a 
distribution of tensile forces in the chains. We verify this 
assumption in section by comparison to the measured 
phonon mode energy in the chains, which reflects the ac- 
tual tensile force in the monatomic wire. All experiments 
to this point were done on gold chains and this is com- 
pared to results for platinum in section IVII which can 
only be qualitative for lack of data available in the liter- 
ature. We conclude (section I V 11(1 that we have a fairly 
good description of Joule heating in these atomic-sized 
conductors and we obtain a realistic estimate for the lat- 
tice temperature at finite bias. 



II. EXPERIMENTAL TECHNIQUES 

The experiments have been performed using the 
Mechanically Controllable Break Junction (MCBJ) 
technique^. By this technique one breaks a wire of the 
desired metal in a controlled fashion with the use of a 
piezo-electric element. The wire is only broken once the 
sample is under vacuum and the temperature in the vac- 
uum chamber is stable at 4.2 K. Thermal contact to the 
sample is provided by a copper finger in order to avoid 
using thermal contact gas. The advantages of the MCBJ 
technique are, firstly, a large mechanical stability due to 
the small mechanical loop connecting the two wire ends. 
This can be of crucial importance when one tries to de- 
termine the stability of monatomic chains. Secondly, the 
cryogenic vacuum without use of thermal exchange gases 
ascertains the absence of gas molecules in the environ- 
ment that could lead to additional cooling of the contact 
and contamination of the metal nanowire. Our relying on 
cryogenic vacuum has the disadvantage that we cannot 
increase the temperature of the experiment much since 
this would lead to a deterioration of the vacuum. The 
gold and platinum wires used here have a purity of at 
least 99.998%. 

Chains of atoms were produced in the following way. 
During the breaking of the wire by use of the piezo- 



element the wire conductance decreases in stepsiS, which 
have been shown to be due to atomic reconfigurationsii. 
The last plateau in the conductance, observed just be- 
fore the conductance suddenly drops to a negligible value 
and the wire is broken, is generally due to a contact with 
the cross section of a single atom.^^. When the length 
of this plateau exceeds the typical length for a single 
atom it can be considered as a measure for the length of 
an atomic chain being formed^ii^. In order to fabricate 
atomic chains of a given length L the process is repeated 
as long as the contact breaks before the plateau (and 
chain) reaches the desired length: The two electrodes are 
pushed together again to reestablish a big contact and 
a new atomic contact is pulled out. When after several 
attempts the desired chain-length is reached the piezo 
movement is halted and measurements on the properties 
of the atomic chain are started. During this procedure 
the conductance of the atomic contact is measured using 
a constant bias voltage of the order of 10 mV. 

The effective length L that we will use as a parame- 
ter below refers to the length of the last plateau in the 
conductance, as was done in previous studies^. The 
length measures the stretching of the contact from the 
point where the conductance first drops to the value of 
a single-atom contact. It is widely accepted that the 
starting point corresponds to a contact of one atom. De- 
pending on the detailed evolution of the chain a plateau 
length of L = 0.5 nm for gold would have two or three 
atoms in a row, and L = 1 nm would have four or five 
atoms. Some recent model calculationsi^ii^ suggest that 
the structure of the contact at the starting point of a 
plateau can be viewed as two electrodes ending in a sin- 
gle atom and touching vertex-atom to vertex-atom which 
would make the chains one atom longer, but this remains 
to be verified experimentally. 

III. RATE OF DECAY OF ATOMIC CHAINS 

One would expect the stability of the monatomic chain 
to be determined by the ratio of the energy barrier 
against breaking, W , and the effective temperature. Toff- 
A population, V, of chains would then show an exponen- 
tial decay as a function of time, t, 

Vit)=r{0)-eM-t/T), r = iexp(VK/fcBTeff), (1) 

where / is the frequency by which the system attempts 
to overcome the barrier. In estimates and simulations 
given below we will take / = 3.5 x 10^^ Hz, a typical 
phonon frequency for monatomic gold chain o^i^^ . 

In order to investigate this stability by a simple experi- 
ment the decay of a population of monatomic gold chains 
was followed as a function of time. The moment a chain 
of the desired length is created is taken as t = and the 
time until the chain spontaneously breaks was stored for 
a set of consecutively formed chains. The evolution of 
the number of surviving chains as a function of time, out 
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FIG. 1: The number of monatomic chains remaining intact 
as a function of time out of an initial population of 100. The 
inset shows the same data on a logarithmic time-axis. The 
results were obtained for monatomic gold chains with an ef- 
fective length of 0.4 nm, at a temperature of 4.2 K and at a 
constantly applied bias voltage of 200 mV. 

of an initial population of 100, is presented in figure ^ 
Clearly, the observed dependence is not described by the 
expected exponential decay. The inset in figure ^ shows 
that it actually more closely follows a logarithm. A log- 
arithm can be obtained if we consider for the chains in 
the population a distribution of barriers, a distribution 
of effective temperatures, or both. Note that the char- 
acteristic time for breaking is influenced by the applied 
bias voltage, which is 200 mV in the example of figure Q 
Near zero bias the lifetime of atomic chains at 4.2 K is 
mainly limited by the patience of the experimentalist. 

IV. A DISTRIBUTION OF BOLTZMANN 
FACTORS 

The distribution of lifetimes for given effective length 
L observed in figure ^ may result from a variation in 
the energy barrier W of the chain, which is likely due 
to variation in the tensile force. From the fact that a 
chain will always break when the tensile force exerted is 
large enough, one can conclude that the effective barrier 
height must decrease with force. As shown by the force 
measurements on monatomic chains performed by Rubio 
et ali^ the force varies strongly while pulling to form a 
chain, and will sensitively depend on the atomic configu- 
rations at the connection points to the banks. Assuming 
a dominant lowest order (linear) term the barrier energy 
can be written as 

W^Wo~a\V\~(3\F\, (2) 

where Wq is the equilibrium barrier energy. We have also 
included a linear decrease of the barrier due to electro- 
migration upon application of a bias voltage V. This 



term is motivated by model calculations by Todorov et 
almi^ Using a single-orbital tight-binding model they 
have shown that the barrier against breaking reduces lin- 
early with voltages up to 1 V. More recent calculations 
based on density functional theory by Brandbyge et al.^^ 
have shown that the electromigration force can be highly 
non-linear. However, for many configurations it will be a 
reasonable approximation. 

Apart from a variation in the energy barrier we should 
also consider the local effective temperature at the point 
contact. For vanishing bias voltages the effective tem- 
perature approaches the bath temperature, Tq. Beyond 
a certain threshold voltage it should rise due to a voltage- 
driven term, Ty. Following reference llTl we describe the 
local effective temperature of the atomic contact as: 

T,s = {T^ + T*)'/\ (3) 

An estimate for the dependence of Ty on voltage can 
be obtained from an earlier paper by Todorov^-^. Three 
terms are included in the estimate: (1) Heating by the 
electrons due to creation of phonons in the atomic wire, 
(2) the cooling by the electrons due to absorbtion of 
phonons, and (3) cooling by the thermal transport of en- 
ergy away from the point contact into the metal banks. 
If only the first two terms are included the tempera- 
ture is simply proportional to the voltage, kBTes — eV , 
where it has been assumed that electrons exchange at 
most a single quantum of vibration energy. Note that 
recent experimental work by Zhitenev et al^ on molec- 
ular junctions indicates that several vibrational quanta 
can be deposited per electron. However, this effect can 
be attributed to the relatively long residing times of the 
electrons inside the molecular bridge^'^. Earlier experi- 
mental work on phonon spectroscopy in monatomic gold 
chains^ shows that multiple phonon scattering is negligi- 
bly small for these ballistic wires, in agreement with the 
assumption by Todorov. 

An effective temperature k^Tcs ~ eV would corre- 
spond to thousands of Kelvins at commonly employed 
bias voltages. However, the third term, the thermal con- 
duction to the banks, limits the heating considerably. It 
is more difficult to make precise estimates, but when the 
behavior of the banks is taken to follow the bulk heat 
conduction and heat capacity we obtain^ 

Tv = l^W\ (4) 

where 7 is a proportionality constant given as 7 = 
60 KV^^/^ nm^^/^ and L is the length of the chain. This 
was qualitatively confirmed by an independent approach 
by Chen et alm^. 

In order to investigate experimentally the effect of bias 
voltage heating on the stability of the chains, taking a dis- 
tribution of wire properties into account, we proceeded 
as follows. First chains of the desired length are pro- 
duced following the procedure described above. This is 
done at a bias voltage of 10 mV, a value which is be- 
low the typical electron- phonon excitation energy-. This 



4 



T ' I ' I ' I ' r 




Break voltage (V) 



T ' 1 ' 1 ■ 1 ' r 




Break voltage (V) 

FIG. 2: A distribution of the number of gold atomic chains 
as a function of bias voltage at which they were observed to 
break. The upper panel shows the experimental and simu- 
lated result for 500 chains with L = 0.5 nm long at a temper- 
ature of 4.2 K, while the lower panel shows similar data for 
chains of L =1.0 nm long. 



ensures that the effective temperature in the point con- 
tact will remain close to the bath temperature. Once a 
chain has been formed we wait for 1 s to see whether the 
chain is indeed stable. Next, the voltage is ramped with 
a constant speed of about IV/s until the chain breaks. 
The bias voltage at breaking is stored and it is reset to 
10 mV. The electrodes are then pushed together again 
to repeat the experiment for the next chain. The re- 
sults for such an experiment are represented in figure |21 
The top panel shows a histogram constructed from the 
measured bias voltages at breaking for a set of 500 gold 
atomic chains with L = 0.5 nm. The distribution clearly 
peaks at 1.2 V with a full width at half maximum of 0.6 V. 
The lower panel shows a similar set of experimental re- 
sults for chains oi L = 1.0 nm in length. In this case 
we observe two peaks in the distribution, one at low bias 
representing a large fraction of chains that break early, 
and a second peak at 0.8 V. A possible dependence on 
the voltage ramp speeds was not investigated, but ear- 
lier work found negligible influence of this parameter in 
atomic contacts^. 

We have suggested above that the breaking of the 



atomic chains can be described as a thermally activated 
process. However, for this we need to verify that the 
energy barrier is not significantly modified by thermal 
expansion due to the rising of the local effective temper- 
ature at high bias voltage. In a test we measured two 
sets of 500 atomic chains. With the first set the bias 
voltage was suddenly increased from 15 to 800 mV once 
the chains had been formed. A large number of chains 
was seen to break following a similar logarithmic time 
dependence as displayed in figure With the second set 
we suddenly dropped the voltage from 300 to 50 mV af- 
ter which we observed not a single chain breaking. Note 
that a higher bias should lead to thermal expansion of 
the electrodesSi and thus to a smaller tensile stress than 
at low bias. These two tests show that the influence of 
the thermal expansion on the energy barrier, if present, 
is dominated by the thermally activated breaking. 

In order to compare the results of figure |21 with the 
proposed model for the variation of the energy barrier 
and effective temperature on bias voltage we need realis- 
tic estimates for the constants in equation |21 We adopt 
the values Wn = 0.738 eV and a = 0.14 eV/V obtained 
in reference [13 from a non-equilibrium tight-binding cal- 
culation for a gold atomic chain. The value for (3 in 
equation[21can be estimated from the results by Rubio et 
ali^ They found that gold atomic chains typically break 
at a tensile force of 1.5 nN. At the point of breaking the 
lifetime, r, of a chain has been reduced to the experi- 
mental time scale, which is 0.1 s in the experiments of 
reference^J The local effective temperature at the point 
contact will be given by the bath temperature in view of 
the low bias voltage of 10 mV used in these experiments. 
For T = 0.1 s and Teff — Tq — 4.2 K we obtain from equa- 
tion ^ an average barrier of 0.011 eV at breaking. The 
barrier of 0.738 eV therefore nearly vanishes under the 
influence of a tensile force of 1.5 nN, resulting in a value 
of /3 = 0.49eV/nN. 

We will now see whether the combination of these val- 
ues give a reasonable description of the experimental re- 
sults of figure [3 We assume here that the only missing, 
i.e. not measured, parameter is the tensile force. The av- 
erage tensile force on an atomic chain can be estimated 
from figure 1 in reference Hll to be 0.95 ± 0.15 nN. We 
performed a numerical simulation of these data based 
on equations ^ to ^ and using a normal distribution of 
forces with a mean value of 0.95 nN and a standard de- 
viation of 0.15 nN, which produces the distribution pre- 
sented in figure [3 The peak position and width agree 
very well with those observed in the experiment. Note 
that all parameters have been obtained from independent 
model calculations and experimental data, so that there 
are no adjustable parameters. For an effective length 
L —1 nm (lower panel) the peak in the numerical simula- 
tion is slightly higher than observed. A more significant 
disagreement is in the width of the distribution, which is 
about twice as large in the numerical result, and also the 
peak at zero bias remains unexplained. Nevertheless, the 
qualitatively agreement obtained using literature values 
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FIG. 3: Position of the maximum in the distribution of break- 
ing voltages (figure|5J for monatomic gold chains as a function 
of their effective length. The bullets represent the experi- 
mental results obtained by fitting a Gaussian curve to each 
population formed by 500 chains. The curve is a numerical 
simulation, as described in the text. 



for the parameters suggests that we have at least cap- 
tured some of the physics involved. Thus, this analysis 
suggests that the broad distribution in breaking voltages 
is mainly caused by the large distribution of tensile forces 
in an ensemble of atomic chains. The processes leading to 
the relatively unstable chains breaking at very low volt- 
ages is unclear as of yet. A possible explanation would be 
a large portion of chains quickly breaking when the bias 
voltage reaches the threshold value of the lowest phonon 
mode in the point contact. 

In order to follow the length dependence in more de- 
tail we have constructed histograms similar to the ones 
shown in figure [21 for many chain-lengths. From the his- 
tograms we extract the position of the maximum of the 
distribution at finite bias and plot this as a function of 
length. This dependence is depicted in figureOand shows 
that the average breaking voltage decreases from 1.2 V at 
L = 0.45 nm to just below 0.6 V at L = 1.25 nm. 

With the same values for the model parameters as used 
above we obtained the curve shown in figure |3| In view 
of the fact that there are no free parameters the agree- 
ment between the model and experiment is surprisingly 
good, although the values at greater length demonstrate 
a larger deviation. The reason for the deviation may be in 
the model or the choice of parameters. However, one as- 
pect that we have ignored and that should become impor- 
tant for longer chains is the contribution of the entropy. 
The larger number of bonds present in longer chains will 
reduce the effective barrier—. 



V. TEST OF THE VARIATION IN TENSILE 
FORCES 

In order to test whether the width in the experimen- 
tal distribution at higher bias voltages is the result of a 
variation in the tensile forces one would like to be able 
to measure the forces in the point contact. The MCBJ 
in its original configuration is not suitable to determine 
these, but it is possible to obtain some information in an 
indirect way using Point Contact Spectroscopy (PCS) — . 
By this technique one probes the differential conductance 
dl/dV as a function of voltage using a small modulation 
superimposed on the dc bias ramp. As soon as the excess 
energy, eV, of the electrons shooting through the point 
contact becomes larger than the energy of a phonon mode 
inside the contact a new channel for back-scattering sets 
in resulting in a sudden reduction of the conductance. 
This technique was applied to atomic contacts by Untiedt 
et alMt and to monatomic gold chains by Agrai't et al&, 
who showed that the observed phonon energy decreases 
when the tensile force on the chain is increased. This is 
as intuitively expected: As the tensile force F increases 
the distance between the atoms a in the chain increases 
leading to a weaker bond and a lower stiffness k — dF/da. 
Since the observed longitudinal phonon mode energy is 
proportional to the phonon frequency changes ac- 
cordingly. 

We make use of this effect in the following experiment. 
First we produce a chain of the desired length via the pro- 
cedure described above. Next, a point contact spectrum 
is recorded by measuring the differential conductance as 
a function of bias voltage in the interval between -50 and 
+50 mV, using a lock-in amplifier. Some chains are lost 
during this procedure due to an increase of the effective 
temperature caused by the increased bias voltage. When 
this happens the data are discarded and a new chain is 
formed. After having measured the differential conduc- 
tance the bias voltage is ramped up similar as in the pre- 
vious experiment, at a rate of 1 V/s. The data are stored 
and analyzed after collecting a large set of curves. For 
all point contact spectra that show a clear phonon signal 
the phonon frequency can be compared to the breaking 
voltage. The data are then sorted according to the ob- 
served phonon frequencies and collected in bins of 1 meV 
wide. For the data within each bin the average breaking 
voltage and its standard deviation are evaluated. 

Figure 0] shows the average breaking voltage as a func- 
tion of the observed phonon frequencies, where the error 
bars on each data point represent the standard deviation. 
For a small fraction of about 2% of the spectra more than 
one phonon frequency could be identified, and these were 
not taken into account in figure^ The break voltage in 
the graph clearly rises as a function of phonon energy, 
just as is expected. For a higher phonon frequency the 
force constant n is higher and the chain is expected to 
be less strained, and thus the energy barrier to break- 
ing is higher. The distribution of breaking voltages at 
each phonon energy is still quite broad. From this ob- 
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FIG. 4: Average breaking voltage of gold chains with an 
effective length of 0.9 nm as a function of the observed phonon 
frequency (bullets). The bars indicate the number of data 
points used to obtain the average (right axis). The standard 
deviation in the break voltage for each phonon frequency is 
represented by the size of the error bar. 



servation one can conclude that the tensile force indeed 
plays the expected role on the stability, but that there 
are still other parameters that influence the survival rate 
of atomic chains at high bias voltage. Such parameters 
could be variations in the shape of the leads giving rise to 
differences in thermal conduction and effective tempera- 
tures. Also, various atomic configurations of the leads at 
the connection points to the chain may give rise to dif- 
ferent equilibrium values for the barrier at zero voltage 
and zero strain, Wq. This may also lead to a variation 
in the dependence of the barrier on bias voltage via elec- 
tromigration forces. 



VI. PLATINUM 

Gold is not unique in forming atomic chains and it is 
now known that this property is shared with Pt and Ir^^. 
In contrast to Au, Pt is a metal with both a significant 
density of s and of d-electrons at the Fermi energy. The 
behavior for Pt could therefore deviate from the results 
predicted by the model presented above, since the the- 
ory was largely based on the single-orbital tight-binding 
studies by Todorov et alM. 

The inset of figure [S] shows a histogram of breaking 
voltages for 500 Pt atomic chains, measured by the same 
method as figure [21 for Au. Since the conductance for 
Pt single-atom contacts or monatomic chains is not as 
sharply defined near a quantum unit of conductance as 
is the case for Au the length of the atomic chain is less ob- 
viously defined. This problem was discussed and success- 
fully addressed in reference and we employ the same 
approach for determining the effective chain length. Al- 
though the results in figure |31 are largely similar to those 
obtained for Au they differ in some respects. First, for 



FIG. 5: The main panel shows the evolution of the maximum 
in the distribution of breaking voltages (illustrated for 500 
chains of length L — 1.0 nm in the inset) for Pt atomic chains 
as a function of length. Filled squares show the data points 
that are in the same interval as measured for Au in figure |3| 
open squares represent data points outside this interval. The 
full curve shows the result of the model calculation described 
in the text. 



an effective length of L = 1.0 nm the position of the max- 
imum is now located at 0.17 V, much lower than found 
for Au, while the width of 0.12 V is only slightly smaller. 
Also, the fraction of chains breaking at very low currents 
is significantly smaller. 

Just as for Au one can follow the position of the max- 
imum in the breaking currents as a function of length. 
Such a graph is presented in figure [S] where the filled 
data points represent data in the same length interval 
that was also measured for Au, while the open symbols 
are at smaller lengths. This is to emphasize that the dif- 
ferences in behavior are less dramatic than a comparison 
of the two figures may suggest. The rapid rise towards 
smaller lengths may not be unique for Pt. We would like 
to use again equations|21through^to describe the exper- 
imental results for the breaking voltage as a function of 
length. However, in contrast to Au most parameters in 
the model cannot be found in the literature. Therefore, 
we take a slightly different approach and start from the 
assumption that all chains break as soon as their mean 
lifetime becomes of the order of the characteristic exper- 
imental time-scale. We further assume that the effective 
temperature at breaking is much larger than the bath 
temperature, so that Tea ~ Ty. Combining equations |3 
-0]we obtain 



_ 1 / WF-a\V\ \ 



(5) 



where Wp — Wq — (3{\F\) represents the barrier against 
breaking averaged over a distribution of forces. We set 
T equal to an experimental time scale of 10 ms, but the 
choice affects the results only logarithmically. Taking 
Wp = 0.23 eV, (note, this is Wq reduced by the averaged 
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strain) 7 = 200KV-i/2njn-i/2 ^ ^ 0.39 eV/V, a set 
of values which is within reasonable limits, the resulting 
curve shown in figure describes the experimental data 
quite well. This set of values is not unique as can be seen 
from inspection of equation|5l the values of a, 7 and Wp 
can be scaled by an arbitrary constant. Yet, we can find 
reasonable values for the constants and obtain a good 
qualitative description of the data using the same model 
as for Au above. For L ^ the break voltage assumes a 
value Wp joL while for large L the break voltage decreases 
as but is expected to saturate when eV approaches 
the characteristic phonon energies of the chain. 

VII. CONCLUSIONS 

We have studied monatomic chains for Au and Pt and 
find a systematic dependence of the characteristic chain- 
breaking voltage vs the chain length. Atomic chains sus- 
tain surprisingly large current densities, but the maxi- 
mum current (or voltage) decreases as the chain length 
grows. We have adopted a model proposed by Todorov 
and coworkers and find a satisfactory agreement with 



the observations. Especially for gold the agreement is 
surprisingly good in view of the fact that we have not 
used any adjustable parameter. Yet, there are devia- 
tions that are not captured by the model and that re- 
quire further study. For Pt the general behavior of the 
average breaking voltage as a function of length can also 
be closely described by the model, although the param- 
eters have not been independently verified. The agree- 
ment obtained inspires confidence in equations 13 and ^ 
and we can now estimate the effective (lattice) tempera- 
ture inside the atomic wire. It rises proportional to the 
square root of of the bias voltage for sufficiently large 
bias and for a monatomic gold chain of length L — \ nm 
a.t V — IV it reaches a temperature of 60 K, which is 
well above the bath temperature of 4.2 K. 
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